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Abstract. The formation and motion of lattice defects such as cracks, dislocations, or 
grain bomidaries, occm'S when the lattice configuration loses stability, that is, when an 
eigenvalue of the Hessian of the lattice energy functional becomes negative. When the 
atomistic energy is approximated by a hybrid energy that couples atomistic and continuum 
' models, the accuracy of the approximation can only be guaranteed near deformations where 

both the atomistic energy as well as the hybrid energy are stable. We propose, therefore, 
, that it is essential for the evaluation of the predictive capability of atomistic-to-continuum 

^ coupling methods near instabilities that a theoretical analysis be performed, at least for 

'""5 some representative model problems, that determines whether the hybrid energies remain 

^NJ stable up to the onset of instability of the atomistic energy. 

We formulate a one-dimensional model problem with nearest and next-nearest neighbor 
interactions and use rigorous analysis, asymptotic methods, and numerical experiments 
to obtain such sharp stability estimates for the basic conservative quasicontinuum (QC) 
approximations. Our results show that the consistent quasi-nonlocal QC approximation 
correctly reproduces the stability of the atomistic system, whereas the inconsistent energy- 
based QC approximation incorrectly predicts instability at a significantly reduced applied 
load that we describe by an analytic criterion in terms of the derivatives of the atomistic 
I " I potential. 

m 
> 

^ 1. Introduction 

^ An important application of atomistic-to-continuum coupling methods is the study of the 

^ quasistatic deformation of a crystal in order to model instabilities such as dislocation for- 

mation during nanoindentation, crack growth, or the deformation of grain boundaries [Is] . 
O In each of these applications, the quasistatic evolution provides an accurate approximation 

>• of the crystal deformation until the evolution approaches an unstable configuration. This 

^ occurs, for example, when a dislocation forms or moves or when a crack tip advances. The 

^ crystal will then typically undergo a dynamic process until it reaches a new stable config- 

^ uration. In order to guarantee an accurate approximation of the entire quasistatic crystal 

deformation, up to the formation of an instability, it is crucial that the equilibrium in the 
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atomistic/continuum hybrid method is stable whenever the corresponding atomistic equi- 
hbrium is. The purpose of this work is to investigate whether the quasicontinuum (QC) 
method has this property. In technical terms, this requires sharp estimates on the stability 
constant in the QC approximation. 

The QC method is an atomistic-to-continuum coupling method that models the contin- 
uum region by using an energy density that exactly reproduces the lattice-based energy 



density at uniform strain (the Cauchy-Born rule) [18 ,21,25 . Several variants of the QC ap- 



proximation have been proposed that differ in how the atomistic and continuum regions are 



coupled [3|[8 ,18 ,26 . In this paper, we present sharp stability analyses for the main examples 
of conservative QC approximations as a means to evaluate their relative predictive properties 
for defect formation and motion. Our sharp stability analyses compare the loads for which 
the atomistic energy is stable, that is, those loads where the Hessian of the atomistic energy 
is positive-definite, with the loads for which the QC energies are stable. It has previously 
been suggested and then observed in computational experiments that inconsistency at the 
atomistic-to-continuum interface can reduce the accuracy for computing a critical applied 
load |8||18(|T9||25]. In this paper, we give an analytical method to estimate the error in the 
critical applied load by deriving stability criteria in terms of the derivatives of the atomistic 
potential. 

Although we present our techniques in a precise mathematical format, we believe that 
these techniques can be utilized in a more informal way by computational scientists to quan- 
titatively evaluate the predictive capability of other atomistic-to-continuum or multiphysics 
models as they arise. For example, our quantitative approach has the potential to estimate 
the reduced critical applied load in QC approximations such as the quasi-nonlocal QC ap- 
proximation (QNL), that are consistent for next-nearest interactions but not for longer range 
interactions. Since the longer range interactions are generally weak, such an estimate may 
give an analytical basis to judging that the reduced critical applied load for QNL with finite 
range interactions is within an acceptable error tolerance. 

The accuracy of various QC approximations and other atomistic-to-continuum coupling 
methods is currently being investigated by both computational experiments and numerical 
analysis (g§g§|g|9)|g[g[T5jl6)|20||23||24|. The main issue that has been studied to date in 
the mathematical analyses is the rate of convergence with respect to the smoothness of the 
continuum solution (however, see [2|[7|[23] for analyses of the error of the QC solutions with 
respect to the atomistic solution, possibly containing defects). Some error estimates have 
been obtained that give theoretical justification for the accuracy of a QC approximation for 
all loads up to the critical atomistic load where the atomistic model loses stability [5 ,23 , but 
other error estimates that have been presented do not hold near the atomistic limit loads. It 
is important to understand whether the break-down of these error estimates is an artifact of 
the analysis, or whether the particular QC approximation actually does incorrectly predict 
an instability before the applied load has reached the correct limit load of the atomistic 
model. 

Two key ingredients in any approximation error analysis are the consistency and stability 
of the approximation scheme. For energy minimization problems, consistency means that 
the truncation error for the equilibrium equations is small in a suitably chosen norm, and 
stability is usually understood as the positivity of the Hessian of the functional. For the 
highly non-convex problems we consider here, stability must necessarily be a local property: 
The configuration space can be divided into stable and unstable regions, and the question we 
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ask is whether the stabihty regions of different QC approximations approximate the stabihty 
region of the full atomistic model in a way that can be controlled in the setup of the method 
(for example, by a judicious choice of the atomistic region). 

In this work, we initiate such a systematic study of the stability of QC approximations. 
In the present paper, we investigate conservative QC approximations, that is, QC approx- 
imations which are formulated in terms of the minimization of an energy functional. In 
a companion paper [6], we study the stability of a force-based approach to atomistic-to- 
continuum coupling that is nonconservative. 

In computational experiments, one often studies the evolution of a system under incre- 
mental loading. There, the critical load at which the system "jumps" from one energy well 
to another is often the goal of the computation. Thus, we will also study the effect of the 
"stability error" on the error in the critical load. 

We will formulate a simple model problem, a one dimensional periodic atomistic chain 
with pairwise next-nearest neighbour interactions of Lennard- Jones type, for which we can 
analyze the issues layed out in the previous paragraphs. It is well known that the uniform 
configuration is stable only up to a critical value of the tensile strain (fracture). We use 
analytic, asymptotic, and numerical approaches to obtain sharp results for the stability of 
different QC approximations when applied to this simple model. 

In Section [2| we describe the model and the various QC approximations that we will an- 
alyze. In Section |4| we study the stability of the atomistic model as well as two consistent 
QC approximations: the local QC approximation (QCL) and the quasi-nonlocal QC approx- 
imation (QNL). We prove that the critical applied strains for both of these approximations 
are equal to the critical applied strain for the atomistic model, up to second-order in the 
atomistic spacing. 

A similar analysis for the inconsistent QCE approximation is more difficult because the 
uniform configuration is not an equilibrium. Thus, in Section [5} we construct a first-order 
correction of the uniform configuration to approximate an equilibrium configuration, and we 
study the positive-definiteness of the Hessian for the linearization about this configuration. 
We explicitly construct a test function with strain concentrated in the atomistic-continuum 
interface that is unstable for applied strains bounded well away from the atomistic critical 
applied strain. 

In Section |6| we analyze the accuracy in predicting the critical strain for onset of instability. 
For the QCL and QNL approximations, this involves comparing the effect of the difference 
between their modified stability criteria and that of the atomistic model. For QCE, since 
the solution to the nonlinear equilibrium equations are non-trivial, we provide computational 
results in addition to an analysis of the critical QCE strain predicted by the approximations 
derived in Section [H 



2. The atomistic and quasicontinuum models 

2.1. The atomistic model problem. Suppose that the infinite lattice eL is deformed 
uniformly into the lattice yp := FeZ, where F > is the macroscopic deformation gradient 
and where £ > scales the reference atomic spacing, that is. 



{yF)i '■= Fie for — oo < ^ < oo. 
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Figure 1. Lennard- Jones type interaction potential. The bond length r^, is 
the turning point between the convex and concave regions of 0. 

We admit 2A^-periodic perturbations u = {u()i(zi from the uniformly deformed lattice yp- 
More precisely, for fixed G N, we admit deformations y from the space 

yp ■■= {y (^^^ ■■ y = yp + u, u ^U], 

where U is the space of 2A^-periodic displacements with zero mean, 

W := {m G : Ue+2N = Ui for i E Z, and J2eL-N+i''^^ ~ O}- 

We set e = 1/N throughout so that the reference length of the periodic domain is fixed. 
Even though the energies and forces we will introduce are well-defined for all 2A^-periodic 
displacements, we require that they have zero mean in order to obtain locally unique solutions 
to the equilibrium equations. These zero mean constraints are an artifact of our periodic 
boundary conditions and are similarly used in the analysis of continuum problems with 
periodic boundary conditions. 

We assume that the stored energy per period of a deformation y G yp is given by a 
next-nearest neighbour pair interaction model, 

N 

£M--=e {<Piyi) + <Piy'i + yi+i)), (i) 

e=-N+i 

where v'^ is the backward difference 

v'^ ■= e-\vi - ve^i) for v eR^, £ e Z, 

and where is a Lennard- Jones type interaction potential satisfying (see also Figure [T| 

(i) 0gC4((O,+oo);M), 

(ii) there exists > such that is convex in (0,r^,) and concave in {r^, +oo), and 

(iii) (f)^^\r) — )■ rapidly as r ^ oo, for A; = 0, . . . , 4. 

We have used the scaled interaction potential, e(f){r/e), in the definition of the stored energy, 
Sa,{y), to obtain a continuum limit as e — )■ 0. Assumptions (i) and (ii) are used throughout 
our analysis, while assumption (iii) serves primarily to motivate that next-nearest neighbour 
interaction terms are typically dominated by nearest-neighbour terms. Note, however, that 
even with assumption (iii), the relative size of next-nearest and nearest neighbour interactions 
is comparable when strains approach r^,. 
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We denote the first variation of the energy functional, £'g^{y)[u], at a deformation y G yp 

by 

^ f)f ( ) ^ 

KiyM:= E ^p^^ = ^ E {4>'iyiH + <P'iy'^ + yi+i)i^'e + ^'i+i)}^ 

£=-N+l l=-N+l 

for u E U. In the absence of external forces, the uniformly deformed lattice y = yp is an 
equilibrium of the atomistic energy under perturbations from U, that is, 

£'^{yp)[u] = forallweW. (2) 

We identify the stability of yp with linear stability under perturbations from the space U. 
To make this precise, we denote the second variation of the energy functional, 
evaluated at a deformation y G 3^^, by 

(•"( 1 d'^£^{y) 
S^{y)[u,v\:= ^ U(Vm 

' ^ (3) 

N ^ ' 

= e U"(y>W£ + <P"iyi + y'i+i)K + K + < 



t=-N+l 

N 



for all u,v eU. The matrix ( 1-%^ ) is the Hessian for the energy functional. We 

say that the equilibrium yp is stable for the atomistic model if this Hessian, evaluated at 
y = yp, is positive definite on the subspace U of zero mean displacements, or equivalently, if 

S"{yp) [u, m] > for all ueU\ {0}. (4) 

In Section [3} Definition [T| we extend this definition of stability to the various QC approxi- 
mation and their equilibria. 

Note that if y = yp, then y[ = F and y[ + y^j^^ = 2F for all i. Therefore, upon defining 
the quantities 

0';:=0"(F), := 0"(2F), and Ap = ^'^ + icP'^^, 

we can rewrite ^ as follows 

N 

S':{yp)[u,u]=e {<I^FK\' + <p2FWi + K+i\'} for ueU. (5) 

e=-N+i 

(We will use Ap later.) The quantities (pp and (f)2p will play a prominent role in the analysis 
of the stability of the atomistic model and its QC approximations and describe the strength 
of the nearest neighbor and next-nearest neighbor interactions, respectively. We similarly 
define the quantities (pQ^ for all A; G N and for all G > 0. For most realistic interaction 
potentials the second-nearest neighbour coefficient is non-positive, — 0; except in the 
case of extreme compression (see Figure [T]). Therefore, in order to avoid having to distinguish 
several cases, we will assume throughout our analysis that F > r,,/2. In this case, property 
(ii) of the interaction potential shows that < 0. 

We also note that, for u G W, both u' and u" are understood as 2A^-periodic chains, that 
is, u',u" G U, where the centered second difference u" eU is defined by 

u" := e~^(M£+i - 2ui + for u eR^, i E Z. 
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For u,v G we also define the weighted 

N 




1 < J9 < OO, 

p = oo, 



as well as the weighted £^-inner product 

N 

{u, v) = e ^ UiVe. 
e=-N+i 

2.2. The local QC approximation (QCL). Before we introduce different flavors of QC 
approximations, we note that we can rewrite the atomistic energy as a sum over the contri- 
butions from each atom, 

N 

£^{y) = e ^ E^{y) where 
e=-N+i 

If y is "smooth," i.e., y[ varies slowly, then E^{y) ^ E^{y) where 

EKy) ■■= '2 Wi) + ^(y'i+i) + ^PC^y'i) + ^iH+i)] = \ [My'i) + My't+i)] , 

and where (pcbi'f') '■= 4>{i^) + 0(2r) is the so-called Cauchy-Born stored energy density. In this 
case, we may expect that the atomistic model is accurately represented by the local QC (or 
continuum) model 

N N 

8^,,{y):=e E^iy) = e My'i)- (6) 

e=-N+i e=~N+i 
The main feature of this continuum model is that the next-nearest neighbour interactions 
have been replaced by nearest neighbour interactions, thus yielding a model with more 
locality. Such a model can subsequently be coarse-grained (i.e., degrees of freedom are 
removed) which yields efficient numerical methods. 

2.3. The energy-based QC approximation (QCE). If y'^ is "smooth" in the majority 
of the computational domain, but not in a small neighbourhood, say, {—K, . . . , K}, where 
K > 1, then we can obtain sufficient accuracy and efficiency by coupling the atomistic 
model to the local QC model by simply choosing energy contributions E^ in the atomistic 
region A = {—K, . . . , K} and E^ in the continuum region C = {—N -t- 1, . . . , N} \ A. This 



approximation of the atomistic energy is often called the energy based QC approximation 21 
and yields the energy functional 

S^M ■.=eY,E'i{y)+eY.E1{y). 

££C eeA 

It is now well-understood [3|j5|[8p5] that the QCE approximation exhibits an inconsistency 
("ghost force") near the interface, which is displayed in the fact that S'^^^lyp) ^ 0. The first 
remedy of this lack of consistency was the ghost force correction scheme [25j which eventually 
led to the derivation of the force-based QC approximation [s] and which we analyze in [6] 
and [T]. 
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2.4. Quasi-nonlocal coupling (QNL). An alternative approach was suggested in [26], 
which requires a modification of the energy at the interface. This idea is best understood in 
terms of interactions rather than energy contributions of individual atoms (see also [8j where 
this has been extended to longer range interactions). The nearest neighbour interactions are 
left unchanged. A next-nearest neighbour interaction (j){e^^{yi^i — yi-i)) is left unchanged 
if at least one of the atoms £ + 1, £ — 1 belong to the atomistic region and is replaced by a 
Cauchy-Born approximation, 

if both atoms belong to the continuum region. This idea leads to the energy functional 

N 

e=-N + l te^qnl teCqnl 

where ^qni = {—K — 1, . . . , K + 1} and Cqni = {— + 1, . . . , A^} \ ^qni are modified atom- 
istic and continuum regions. The QNL approximation is consistent, that is, y = yp is an 
equilibrium of the QNL energy functional. The label QNL comes from the original intuition 
of considering interfacial atoms as quasi-nonlocal, i.e., they interact by different rules with 
atoms in the atomistic and continuum regions. 

3. Stability of Quasicontinuum Approximations: Summary 
In this section, we briefly summarize the our main results. 

We begin by giving a careful deflnition of a notion of stability. Our condition is slightly 
stronger than local minimality, which is the natural concept of stability in statics. However, 
an analysis of local minimality alone is usually not tractable. Moreover, for the deformations 
that we consider, our deflnition is in fact sufficiently general. 

Definition 1 (Stable Equilibrium). Let E : yp — > M U {+00} . We say that y G yp is a 
stable equilibrium of £ if S is twice differentiable at y and the following conditions hold: 

(i) £'{y)[u] = for allu eU, 
(11) £"{y)[u, m] > for allueU\ {0}. 

// only (i) holds, then we call y a critical point of S. 

We focus on the deformation yp = {Fie)e^z and ask for which macroscopic strains F this 
deformation is a stable equilibrium. We know from ^ that yp is a critical point of the 
atomistic energy Sg., and it is easy to see that yp is also a critical point of the QCL energy 
Sqci and of the QNL energy S^^i. Our analysis in Section |4] gives the following conditions 
under which yp is stable: 

1. yp is a stable equilibrium of if and only ii Ap — e^vr^^j^ + O(e^) > 0; 

2. yp is a stable equilibrium of S^ci if and only ii Ap > 0; 

3. yp is a stable equilibrium of S^^i if and only ii Ap > 

where we recall that Ap = (f)p + 402^, is the continuum elastic modulus for the Cauchy- 
Born stored energy function (pcbij) = 0('^) + 0(2r). Points 1., 2., and 3. are established, 
respectively, in Propositions |2} [3} and |4j 

If we envision a quasistatic process in which F is slowly increased, then we may wish to 
find the critical strain F* at which yp is no longer a stable equilibrium (fracture instability). 
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If we denote the critical strains in the atomistic, QCL, and QNL models, respectively, by 
F*, and F*^y, then 1.-3. imply that (cf. Section g 

\F: - F*J = 0{e') and \F: - F*J = Oie'). 



For the QCE approximation defined in Section 2.3| the situation is more complicated. The 
occurrence of a "ghost force" in the QCE model implies that yp is not a critical point of S^ce, 
and consequently, we will need to analyze the stability of the second variation S'^^^{yqce,F) 
where yqce,F 7^ yp is an appropriately chosen equilibrium of S^ce- Since yqce,F solves a 
nonlinear equation, we will replace it by an approximate equilibrium in our analysis in 
Section [5] where we obtain the following (simplified) result: 

4. For yqce,F to be a stable equilibrium of Sqce it is necessary that 

20'; +2|0';|2+'^^'^>^' 

where S = max {|0(^)(2F)/0"(F)| : j = 1,2,3} is assumed to be small. 

We remark that 4. gives only a necessary but not a sufficient condition for stability of the 
QCE equilibrium yqce,F, which, moreover, depend on assumptions on the parameter 6. We 
refer to Remark |3] for a careful discussion of the role of S. 

If we let F*^g denote the critical strain at which 4. fails (ignoring the 0{6^) term), then 
we obtain 

\F:-F*^J = Oil), 

which suggests that the QCE method is unable to predict the onset of fracture instability 
accurately. In Section |6| we confirm this asymptotic prediction with numerical experiments. 

We have shown in [17] that the stability properties of the ghost force correction scheme 
(GFC) can be understood for uniaxial tensile loading by considering the stability of the QC 
energy 

^gfcj/f (y) := ^qcc(y) - ^qce(l/F)(z/ " Z/f) for all y E yp- (7) 

We note that S'gfc, ypiVF) = 0, so yp is an equilibrium of the £gic,yF energy under perturbations 
from lA. We can therefore analyze the stability of Sgicypiy) at yp by studying the Hessian 
^qcehjp) = ^gicypiyp)- We show in Remark[2]in Section 



2.3 



that 



5. yp is a stable equilibrium of £gic,yp if and only ii Ap + Xk<P2f > where ^ < Xr < 1. 

We give analytic and computation results in Sections [5] and |6] showing that the ghost force 
correction scheme can be expected to improve the accuracy of the computation of the critical 
strain by the QCE method, that is. 



qcc qcc ^ a ' 



where F^^.^ is the critical strain at which S'^^^{yp) = £gic,yphjF) is no longer positive definite, 
but the error in computing the critical strain by the GFC scheme is still 



\F:-Fy^\ = o{i) 
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4. Sharp Stability Analysis of Consistent QC Approximations 

In this section, we analyze the stabihty of the atomistic model and two consistent QC 
approximations: the local QC approximation and the quasi-nonlocal QC approximation. In 
each case, we will give precise conditions on F under which yp is stable in the respective 
approximation. The inconsistent energy-based QC approximation (QCE) is analyzed in 
Section [5j The corresponding result for QCE is less exact than for QCL and QNL, but 
shows that there is a much more significant loss of stability. 

4.1. Atomistic model. Recalling the representation of S'J^^yp) from ([s]) and noting that 



|/i / |2 Ol'|2_i_o|' |2 \ ' I \ 2 
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we obtain 

N N 



'-\u'.-u:\') 



e=-N+l E=-N+l 

N N 

i=-N+i e=-N+i 
= Ap\\u% + i-e'<P'ip)\\u%. (9) 

To quantify the influence of the strain gradient term, we define 

jjLr := ml - — — . 
^ew\{o} ll^'lls 

Since u is periodic, it follows that u' has zero mean. In this case, the eigenvalue /i^ is known 



to be attained by the eigenfunction ■?/'^ = sin(£:£7r) and is given by 27, Exercise 13.9] 



= (10) 

Since sin(t) = t + 0(1?) as t \ 0, it follows that /i^ = tt + O(e^) as e \ 0. Thus, we obtain 
the following stability result for the atomistic model. 

Proposition 2. Suppose (p'J^p < 0. Then yp is stable in the atomistic model if and only if 
Ap — e'^fil(f)2p > 0, where fi^ is the eigenvalue defined in (10). 



Proof. By the definition of /ig, and using (|9]), we have 



inf S'^{yF)[u,u] = Ap-e cji'^p inf \\u"\\p = Ap - e fi^cp'^p. □ 



||m'||^2=i II«'I|£2=i 



4.2. The Local QC approximation. The equilibrium system, in variational form, for the 
QCL approximation is 



TV 



£UyM=^ E (0'(z/;) + 20'(2y;))w^ = O forallweW. 



t=-N+l 
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Since u' has zero mean, it follows that y = yp is a critical point of Sqci for all F. The second 
variation of the local QC energy, evaluated at y = yp, is given by 

N 

e=-N+i 



Thus, recalling our definition of stability from Section 2.1, we obtain the following result 



Proposition 3. The deformation yp is a stable equilibrium of the local QC approximation 
if and only if Ap > 0. 

Comparing Proposition [3] with Proposition [2] we see a first discrepancy, albeit small, be- 
tween the stability of the full atomistic model and the local QC approximation (or the 
Cauchy-Born approximation). In Section |6] we will show that this leads to a negligible error 
in the computed critical load. 

4.3. Quasi-nonlocal coupling. By the construction of the QNL coupling rule at the in- 
terface, the deformation y = yp is an equilibrium of Sqri\ |26j. The second variation of Sqni 
evaluated at y = yp is given by 

N 

i=-N+i eeA^ni 

+ e J2 402V(iKI' + ^l%+il')- 

We use ([s]) to rewrite the second group on the right-hand side (the nonlocal interactions) 
in the form 

K+l K+l 

e=-K-i e=-K-i 

to obtain 

N K+l 

^qni{yF)[u,u] = e Y Apluf + e (-^^^2f)K\'^- 
e=-N+i e=-K-i 

Except in the case K G {A^ — 1, A^}, it now follows immediately that yp is stable in the QNL 
approximation if and only ii Ap > 0. 

Proposition 4. Suppose that K < N — 1 and that (p^p < 0, then yp is stable in the QNL 
approximation if and only if Ap > 0. 

5. Stability Analysis of the Energy-based QC approximation 
We will explain in Remark |2] that there exists | < A^^- < 1 such that 

^qceiVF) is positive definite if and only if Ap + Xk'P'2p > 0- 
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However, yp is not a critical point of £qcc, so we must be careful in extending the previous 
definition of stability to the QCE approximation. We cannot simply consider the positive- 
definiteness of E'^^^iyp). Instead, we analyze the second variation S'^ceiyqce,F) where yqce,F ^ 
yp solves the QCE equilibrium equations 

d£ 

^(l/qce,F)=0 for£=-iV+l,...,iV, 

or equivalently 

^;ce(2/qce,F)M = for all u G W. (11) 

We will see that, when the second- neighbour interactions are small compared with the first 
neighbour interactions (which we make precise in Lemmajs]), there is a locally unique solution 
2/qce,F of the equilibrium equations, which is the correct QCE counterpart of yp. We will 
then derive a stability criterion for the equilibrium deformation yqce,F- 

In Proposition 6] below, we derive an upper bound for the coercivity of £'^^^{yqce,F)[u, u] 
with respect to the norm ||^^'||^2. Even though the derivation of this upper bound is only 
rigorous for strains bounded away from the atomistic critical strain, it clearly identifies a 
source of instability that cannot be found by analyzing, for example, £'^^^{yp). Moreover, 
we will present numerical experiments in Section [6] showing that the critical strain predicted 
in our following analysis gives a remarkably accurate approximation to actual QCE critical 
strain. 

In Section [6} we consider the critical strain for each approximation, namely the point at 
which the appropriate equilibrium deformation (either yp or yqce,F) becomes unstable. We 
will see later in this section, as well as in Section [6] that predicting the loss of stability for 
the QCE approximation using yp greatly underestimates the error in approximating the 
atomistic critical strain by the QCE critical strain. 

Due to the nonlinearity and nonlocality of the interaction law, we cannot compute yqce,F 
explicitly. Instead, we will construct an approximation ^qce,F which is accurate whenever 
second-neighbour terms are dominated by first-neighbour terms. In the following paragraphs, 
we first present a semi-heuristic construction, motivated by the analysis in [4], and then a 
rigorous approximation result, the proof of which is given in Appendix [B] 



In (26) in the appendix, we provide an explicit representation of S'^^^^. Inserting y = yp, 
we obtam a variational representation of the atomistic-to-continuum interfacial truncation 
error terms that are often dubbed "ghost forces," 

:= —(t>'2p{g\ u') for u eU, 



where 



-|, i= -K -1,K + 2, 
9e={ i i=-K + l,K, (13) 
0, otherwise. 



We note that (12) makes our claim precise that yp is not a critical point of Sqce- 

Motivated by property (iii) of the interaction potential (p, we will assume that the param- 
eters 

and S,_:^=f^ 
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are small, and construct an approximation for ?/qcc,F which is asymptotically of second order 
as ^1,52 — !■ 0. Although such an approximation will not be valid near the critical strain for 
the QCE approximation, it will give us a rough impression how the inconsistency affects the 
stability of the system. We note that 6i is scale invariant since we used a scaled interaction 
potential, e(j){r/e)^ in our definition of the stored energy ([T]). 

A non-dimensionalization of (12) shows that Uqce.F = Vf + 0{6i). If 6i is small, then we 
can linearize (11) about yp and find the first-order correction yn^ G yp, which is given by 



^qceiyp) [yim " u] = S' (yp) [u] = (^'^f^q', for all u^U. 



(14) 



We note that this linear system is precisely the one analyzed in detail in [3] . However, instead 
of using the implicit representation of yun — yp obtained there, we use the assumption that 
^2 is small to simplify ( 14 ) further and obtain a more explicit approximation. 

Writing out the bilinear form E'^^Jyyp^\ii, n] explicitly (using (28) as a starting point) gives 

N A'-l 



+ 102V 



i=0 
+ 102V 



N 

E 

•.=K+2 



4:(f)''Ju 



I |2 



(15) 



u 



K+l 



U 



K+2\ 



+ f402V|w'i 



2Fl"X+ll 



where we have only displayed the terms in the right half of the domain and indicated the 
terms in the left half by dots. Ignoring all terms involving (p'/^p, which are of order 82 relative 
to the remaining terms, we arrive at the following approximation of (14): 

0F((yqce,F " Z/f)', u') = (^^F^Q ^ all U G W, 

the solution of which is given by 

The following lemma makes this approximation rigorous. A complete proof is given in 
Appendix [B| 



Lemma 5. If 8\ and 62 are sufficiently small, then there exists a (locally unique) solution 

yqce,F of (11) SUCh that 

\\{yqce,F - yqce,F)'||<?°° < C{6l + 51^2), 

where C may depend on cf) (and its derivatives) and on F, hut is independent of e. 

From now on, we will also assume that ^3 := 02^/0'^ is small, and combine the three small 
parameters into a single parameter 

5 := max(|5i|, l^al, 

We will neglect all terms which are of order 0(5^). A careful discussion of the parameter 5 
and the validity of the asymptotic analysis is given in Remark [3j 

In the following, we will again only show terms appearing on the right half of the domain. 
Our goal in the remainder of this section is to obtain an estimate for the smallest eigenvalue 
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of £^qj,g(^qce,F)- Using ([28|, we can represent S'^^^{yqce,F) as 

K-2 N 



e=0 i=K+3 

+e{<p"{F + ISi) + 3<P"{2F + i^i)} \u'k\^ 
+e{(P'^ + 0"(2F - i^i) + 0"(2F + i^i) + 20'2V} 
+e{(P"{F - + 0"(2F - i^i) + 40"(2F - 6^)} \u'k^,\ 



We expand all terms containing 6i and neglect all terms which are of order 0(5^) relative 
to (pp, which is the order of magnitude of the coefficient of the diagonal term of £^^.g{yc^ce,F) ■ 
For example, we have, for some G (0, 1), 

~ = ^ + jr> (2^1) = ^ + 0{d3di), 



y^F H^F 



as 61,63 — 7- 0. Thus, the 0{6i) perturbation of a second-neighbour term will not affect our 
final result. On the other hand, expanding a nearest neighbour term gives 



nF+j6^_ 



^ + §i¥i) + o{6!) 



YF 'TF 

as (5i — )■ 0. Proceeding in the same fashion for the remaining terms, we arrive at 

K-l K~l N 

£=0 i=0 e=K+3 

+ e{AF + il6i(pF - 02f)}KI' + (16) 

+ e{AF - - 02F)}|WW2P - ^^l<p2F{\u'k\'^ + K'+lT} 

+ O(0';5ln'|||). 

Clearly, our focus must be the coefficients of the terms jw'^^P and |M'^_,_2p) and in particular, 
on the quantity 

IXA"' A" _ 4>"f4>'2F ~ '^4'"f4>'2F /^^7^ 
i^OlCpp - (p2F - 1^1 • \^') 

Depending on the sign of \6i(t)"p — < 0, we see that the "weakest bonds" are either 
between atoms K — 1 and K (as well as —K + 1 and —K) or between atoms K + 1 and 
+ 2 (as well as -fsT - 1 and -K - 2). 

If ^6i(j)'p — < 0) insert the test function w gU, defined by 



[y^y/^ i = K, 

0, otherwise. 
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into (16) to obtain 



ll«'ll^2=l 



Af\1 + 



2F 



+ 



;i8) 



Note that the constant 2 in front of 4'p4>2F was replaced by 5 due to the strain gradient terms 
in (16) which shghtly stabihze the system. 

If 2^i(f)p — (j)'2F > 0, we use the alternative test function w gU, defined by 



(i£-i)i/2, e = K + 2, 

0, otherwise, 



(19) 



to test (16), which gives 



lh'IL2=i 



A, 



^'i^<P'2F 



%<^'Lf 



2Ap<^% 



+ 



0(5^)}. 



(20) 



In this case, only a single strain gradient term affects the final result, and therefore this 
correction is only small. 

Due to the stabilizing effect of the strain gradient terms for our perturbation, the right 



hand sides of (18) and (20) might both be bounded below by Ap^ so our estimate will involve 



a min over three terms. Recalling that ?/qce,F = yqce,F + C'((5^), we obtain the following result: 



Proposition 6. There exist constants 6 and C, which may depend on (j) and its derivatives 
and on F but not on e, such that, if S < 6, then 



inf ^qce(2/qce,F) [u, u] < <j)"p (mini 1 + 



Il«'llf2=l 



± 



2|0';|2 



3(j)'^p\ Ap , ^ 2 

Iir]^^\+CS j. (21) 



Proof. The bounds ( 18 ) and (20 ) are rigorous provided 6 is sufficiently small so that F — ^6i is 



bounded away from zero. Moreover, if 6 is sufficiently small, then Lemma |5] gives a rigorous 
bound for the error ||(2/qce,F ~ yqce)'||^°° which only adds an additional 0(5^) error to the 
estimate. □ 



For typical interaction potentials, we would expect that <; 
that 02F > 0) ^6 have already postulated that (pp > 



< 
and 



'as (l)p is decreasing), 
< 0. Thus, the 



two terms in the numerator of the right hand side of (17) have opposing sign and may, in 



principle even cancel each other. However, we have found in numerical tests that for typical 
potentials such as the Morse or Lennard- Jones potentials the first term is dominant, that 



IS 



' 2 



- i02F < 02^ and 



mm<! 1 H -^r- ± 



Af 



1 + 



^2F ^ m'2F 



^'f 
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in Proposition [6j 

Remark 1. Proposition |6] as well as the subsequent discussion clearly shows that the 
spurious QCE instability is due to a combination of the effect of the "ghost force" error and 
of the anharmonicity of the atomistic potential. □ 

Remark 2. A variant of the analysis presented above shows that S'^^^lyp) is positive definite 
if and only if Ap + Xk4>2f > where ^ < A/^ < 1. The lower bound can be obtained using 
the test function (19) in the bilinear form S'^^^{yp)[u, u] given explicitly by (15), while the 



upper bound can be obtained from the estimate 

^qcciyF)[u, U] > i^F + 02f)II'"'II| ^11 M G W, 

which also follows from (15) (see also [Hj Lemma 2.1]). Thus, the lower bound is related 



to the second term in (21) which we have noted above is generally greater than the first 
term, and we can conclude that the critical strain for QCE obtained by linearizing about yp, 
rather than the equilibrium solution yqce,F, significantly underestimates the loss of stability 
(see also Figure [2|. 



The study of the positive-definiteness of S''{yp) is relevant to the stability of the ghost- 



force correction iteration and is discussed in more detail in 17 . □ 



Remark 3. While our rigorous results. Lemma [5] and Proposition [6} are proven only for 
sufficiently small 5, one usually expects that such asymptotic expansions have a wider range 
of validity than that predicted by the analysis. For this reason, we have neglected to give 
more explicit bounds on how small 5 needs to be. 

Nevertheless, a relatively simple asymptotic analysis such as the one we have presented 
cannot usually give complete information near the onset of instability. Our aim was mainly 
to demonstrate that the inconsistency at the interface leads to a decreased stability of the 
QCE approximation when compared to the full atomistic model or the consistent QC ap- 



proximations. We will see in Section |6j that, if we use (21) to predict the onset of instability 
for QCE, then we observe a fairly significant loss of stability of the QCE approximation when 
compared to the full atomistic model. In numerical experiments, we will also see that the 



prediction given by (21) is qualitatively fairly accurate for the Morse potential for a range 



of parameters that explores the dependence of our results on 5. □ 

6. Prediction of the Limit Strain for Fracture Instability 

The deformation yp G yp is an equilibrium of the atomistic energy for all F > 0. However, 
it is established in Proposition [2] that yp is stable if and only if F < F* where F* is the 
solution of the equation 

MF:) := 0"(F:) + (4 - eV?)0"(2F:) = 0. (22) 

We call F* the critical strain for the atomistic model. The goal of the present section is to use 
the stability analyses of the different QC approximations in Sections |4] and |5] to investigate 
how well the critical strains for the different QC approximations approximate that of the 
atomistic model. 

In order to test our predictions against numerical values, we will use the Morse potential 
0^(r) = e-^^C'--!) - 2e-"('^-i) = (e'^^"-^) - 1)^ - 1, (23) 
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02 03 04 05 06 07 


0ij 


Cerr(0) 


1.0877 0.3796 0.1339 0.0485 0.0177 0.0065 


0.0635 



Table 1. Numerical values of the error constant Ccrr(0) defined in (24), for 
various choices of 0. 
where a > 1 is a fixed parameter, and the Lennard- Jones potential 

6.1. Limit strain for the QCL and QNL approximations. The critical strain F* for 
the local QC approximation as well as the QNL approximation (cf . Propositions |3] and |4]) is 
the solution to the equation 

MK) :=0"(f;) + 40"(2f;) = o. 

We note that the critical strain F* for the QCL and QNL models is independent of 
which is convenient for the following analysis. Inserting F* into (22) gives 

Mf:) = Mf:) - eV?0"(2F;) = -eV?0"(2F;), 

and hence 

V'.(F:)-^,(F:)=eVe0"(2F;). 
A linearization of the left-hand side gives 



- f:) = eVe0"(2i^c*) + o{\f: - f. 

Noting that ^^^(^c ) = ^c(^c ) + 0(6^), we find that the relative error satisfies 

+ 0{e') 



f: 


-f: 




vr20"(2F;) 


Fo 


— TP* 

c 




(0'"(F*) + 80'"(2F^*))(Fo - F*) 



(24) 



■=e^CU(t>) + 0{e^), 

where Fq is the energy-minimizing macroscopic deformation gradient which satisfies 



dF 



-(Fo) = 0'(Fo) + 20'(2Fo) = O. 



In Table [T] we display numerical values of Ccrr(0) for the Morse potential = 0a, with 
a = 2, ... ,7, and for the Lennard-Jones potential = 0ij. We observe that the constant 
decays exponentially as the stiffness increases, and that it is moderate even for very soft 
interaction potentials (Ccrr(02) ~ 1.0877). 

6.2. Limit strain for the QCE approximation. In Section|5| we have computed a rough 
estimate for the coercivity constant of the QCE approximation. We argued that, for as long as 
the second neighbour interaction is small in comparison to the nearest neighbour interaction, 
we have the bound 



inf ^qcc (l/qce.F) [u,u]< | 



1 + 



lh'IL2=i 



30|^ 
2 0'; 



+ 



0;V2F 

2|0';P 



Even though this bound will, in all likelihood, become invalid near the critical strain, it is 
nevertheless reasonable to expect that solving 



+ i02F + 



^F02F 

20'^ 



0, 



(25) 
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Figure 2. Critical strains F, 



qce' qce' 



qce' 



F* and the equilibrium strain 
Fq, computed for the Morse potential (23) with varying a. The critical strains 
for the QCE Hessian, F*^^, are computed with = 40 and K = 10. The 
approximation, F*^^, is computed using the asymptotic approximation (25). 
The strain F^^^ is the critical strain at which £'^^^{yp 
definite. 



is no longer positive 



will give a good approximation for the exact critical strain, F*^^. The latter is, loosely 
speaking, defined as the maximal strain F > for which a stable "elastic" equilibrium of 
£^qce exists in 3^^;^. A deformation y can be called elastic if = 0(1) for all £, as opposed to 
fractured if = 0{N) for some Eq. 

We could use the same argument as in the previous subsection to obtain a representation 
of the error; however, since F*^^ depends only on F but not on e we can simply solve for F*^^ 
directly. 



qce 



For the Morse potential (23), with stiffness parameter 2 < a < 7, we have computed both 
(for = 40, i^' = 10 as well as for A^ = 100, K = 20) and F*^^ numerically and have 
plotted these critical strains in Figure [2| comparing them against Fq and F* . We have also 
included the critical strain F^^^, below which S'^^^^yp) is positive definite, to demonstrate that 
it bears no relation to the stability or instability of the QCE approximation. We discuss F^^^ 
in detail in jl^ where we argue that it describes the stability of the ghost-force correction 
scheme. 

In Figure [3} we plot the relative errors 



a I— J- 



qce 



'a)-F:(a] 



F*(a)-Fo(a) 



and a h-> 



qce 



'a)-F*(a) 



F*(a)-Fo(a) 



We observe that the prediction for the critical strain, as well as the prediction for the 
relative error, obtained from our asymptotic analysis is insufficient for very soft potentials but 
becomes fairly accurate with increasing stiffness. In particular, it provides a good prediction 
of the relative errors for the critical strains for a > 3.5. 

For a correct interpretation of our results, we must first of all note that the relative errors 
for the critical strains decay exponentially with increasing stiffness a. While, for small a 
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Relative Errors of Critical Strains 



25 -< 



g 
m 



0.25 



0.04 













^F^,N = 100, K = 20 




qce 




qce 





a 



7 



Figure 3. Relative errors of the critical strains (computed and predicted) 
for the QCE approximation against the critical strains of the QCL/QNL ap- 
proximation. The errors are computed explicitly for = 40, K = 10 as well 
as for = 100, K = 20, using the Morse potential (23) with varying a. These 
two curves are very close and may be hard to distinguish. Additionally, we 
show the critical strain for loss of positive definiteness of E'^^^iyp), which does 
not predict the loss of stability that the QCE experiences correctly for any 
parameter value. 



(soft potentials) the error is quite severe, one could argue that it is insignificant (i.e., well 
below 10%) for moderately large a (stiff potentials). However, our point of view is that, by 
a careful choice of the atomistic region one should be able to control this error, as is the 
case for consistent QC approximations such as QNL. For the QCE approximation, this is 
impossible: the error in the critical strain is uncontrolled. 

Conclusion 

We propose sharp stability analysis as a theoretical criterion for evaluating the predic- 
tive capability of atomistic-to-continuum coupling methods. Our results show that a sharp 
stability analysis is as important as a sharp truncation error (consistency) analysis for the 
evaluation of atomistic-to-continuum coupling methods, and provides a new means to dis- 
tinguish the relative merits of the various methods. Our results also provide an approach to 
establish a theoretical basis for the conclusions of the benchmark numerical tests reported 
in particular for the poor performance of the QCE approximation in predicting the 



m 



19 



movement of a dipole of Lomer dislocations under applied shear. 

Of course, the simple one-dimensional situation that we have considered here cannot nearly 
capture the complexity of atomistic-to-continuum coupling methods in 2D/3D. Even the 
much simpler question of whether QCL (the Cauchy-Born continuum model without cou- 
pling) can correctly predict bifurcation points becomes much more difficult since it is possible, 
in general, that the stability region for the Cauchy-Born model is much larger than that for 
atomistic model 10 . However, in many interesting situations this effect does not occur [13] , 



and it is an interesting question to characterize these. Concerning the stability of the coupling 
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mechanism, no rigorous results are available in 2D/3D. Until such an analysis is available, 
we propose that careful numerical experiments should be performed, which experimentally 
investigate the stability properties of atomistic-to-continuum coupling methods. 



Appendix A. Representations of S'^^^ and S'^^^ 

Our aim in this section is to derive useful representations for the first and second variations 
^qceiy) ^iid E'^^^iy) of the QCE energy functional. For notational convenience, we will only 
write out terms in the right half of the domain {— A + 1, . . . , A}, indicating the remaining 
terms (which can be obtained from symmetry considerations) by dots. For example, we write 

N K-l N 

e=0 e=0 i=K+2 

+ IHv'k + v'k+i) + l<t>{y'K+i + y'K+2) + |0(2z/Wi)- 

The first variation is a linear form on W, given by 

AT K-l 

£=0 £=0 

AT 

+ |'/''(2|/Wi)(2«Wi) + ^ E '^'(2l/£)(2«^)- 
Collecting terms related to element strains m^, we obtain 

K-\ 

+ + y'x) + ¥{y'K + yWi)} < 
+ i*/-'!?/^ + y'K^i) + ^'/•'(i/Wi + + '/•'(21/Wi)} <+i (26) 

+£{0'(l/^+2 + \<\>[yK+x + + 20'(2y^+2)} 

AT 

+£ E {0'(y;) + 20'(2y^)}x^i. 
Similarly, the second variation can be written in the form 

N K-l 

+l^"{y'K + y'K+i)K + «Wir + l<l>"iy'K+i + y'K+2)K+i + y^K+2? (27) 

N 
i=K+2 
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Using ([8j) to replace all second- neighbour terms in (27), we obtain the alternative represen- 
tation 

K-1 

2 



7^+21 
N 

.1 |2 



-.^ E + - ¥'{<P'\y'K + + + y'K+2)WK+i?]- 

While somewhat unwieldy at first glance, this representation is particularly useful for the 
stability analysis in Section [5j 

Appendix B. Proof of Lemma [5] 
In this section, we complete the proof of Lemma [5] which was merely hinted at in the main 



text of Section |5j Recall that ^qce,F = yp + ^ig where g is given by (13), and recall, moreover, 
that ^qce solves the linear system 

</>F((yqce,F - ypy.u') = ^'^p{g',u') = -S'^^Mlu] for all ueU. (29) 

Our strategy is to prove that yqce,F has a residual of order 0(5^ + 51^2) and that i^qce(yqce,F) is 
an isomorphism between suitable function spaces. We will then apply a quantitative inverse 
function theorem to prove the existence of a solution yqcc,F of the QCE criticality condition 



(11) which is "close" to yqce,F- Before we embark on this analysis, we make several comments 
and introduce some notation that will be helpful later on. 

To ensure that S^ce is sufficiently differentiable in a neighbourhood of yqcc,F we only need 
to assume that F > and that 61 is sufficiently small, e.g., 61 < F. In that case, S^ce is 
three times differentiable at y for any y e yp such that \\y' — yqce,F\\i°° < ^^i- 

We will interpret S'^^^ as a nonlinear operator from to U~^'°° which are, respectively, 
the spaces lA and U* endowed with the Sobolev-type norms, 

||m||2^i,oo = IIm'II^oo for u eU, and ||T||2^-i,oo = sup T[v] for T G U*. 



Consequently, for y G yp, ^'Leiv) can be understood as a linear operator from to U 



— 1,00 



Our justification for defining yqcc,F as we did in (29) is the bound 

\^qcci.yF)[u,v] - (j)"p{u' ,v')\ < (t)"FCi52\\u'\\i,^\\v'\\ii for all u,v eU, (30) 
where Ci = 5, which follows from (15). We can formulate this bound equivalently as 

W^qceiyp) - 0F^i||l(W1>-, < 4>fCiS2, (31) 

where Li : W — > is given by 

Li{u)[v] = {u',v') for all u,v eU. 
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We also remark that Li : — )• U is an isomorphism, uniformly bounded in N, more 
precisely, 

11-^1 ^\\l{u-'^'°°, u^'°°) ^ 2. (32) 

This result follows, for example, as a special case of |23^, Eq. (36)] or [7| Eq. (5.2)], and is 
also contained in j6j. 

We are now ready to estimate the residual of yqce,F- Expanding S'qceiyqce,F) to first order 
gives 

^qcc(^qcc,F)H = {^qccd/^) H + ^l^qccl?/^) ^] } 

+ -^i / X^qceiyF + thg) [g, v] - Sq^^iyp) [g, v] } dt. 
Jo 



We will estimate the two groups on the right-hand side of (33) separately. Using (12) and 



(30), we obtain 

|^qce(Z/F)H + 5l^qce(l/F)[^,f]| = ^1 | - 0f(^'>^') + ^qco (Z/f) ^^] | 

< 0p'Ci5i52||5''||f°° ll^'ll^i for all v eU. 



(34) 



To estimate the second group in (33) we simply use the regularity of the interaction 
potential (we assumed that G C^(0,+oo)) and Holder's inequality to obtain 



S^,,{yp + t6,g)[g,v]-S'^Jyp)[g,v]\ < <P';.C2t6,\\g'\\UW\Ui 



(35) 



where ((/>'^C2) is a local Lipschitz constant for (p", that is, there exists a universal constant C2 

max(|0'"(F + r)|,|0'"(2(F + r))|) 



C2 = C2 sup 

.1 



\r\<2^l 

In particular, if 5i is sufficiently small then we may assume that 



C2 = 2C2 



max 



I/// 1 I J,/// 



Inserting (35) and (34) into (33), and using the fact that = |, we obtain the 



U ^'°°-residual estimate 

ll^qce(^qcc,F)||w-l>°° < 0f(^Ci(5i(52 + |c2(5i). 

Next, we estimate ||i^qce(yqce,F)"lL(w-i.- 



^^1,00). Using (31) and a similar argument as for 



(35) gives 



l^qce(^qce,F) " 1| L(Wl>- , W"!-) < ll-^qcel^qce.F) " ^qceiV f)\\ L{U^'-° , U-^°°) 

+ W^'lcciVF) - 0F-^l|U(Wi.°°, 
< 0f(^C25i + Ci52). 



Moreover, from (32), we deduce that 



^F 
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A standard result of operator theory states that if X, Y are Banach spaces and T, S : X ^ 
Y are bounded hnear operators with T being invertible and satisfying US' — T|| < 1/||T^^||, 
then S is invertible and 

\\s-'\\ < . 

In our case, setting T = 4>pLi and S = S'^^^{yqce,F) , this translates to 

2 

IIC(w)-1l(..-.oc,^..) < _ - c,5,y 

provided that the denominator is positive. Thus, for 61,62 sufficiently small, we obtain the 
bound 

// - 1 4 

We now apply the following version of the inverse function theorem. 

Lemma 7. Let X,Y be Banach spaces, U an open subset of X , and let F : U ^ Y be 
Frechet differentiable. Suppose that xq G f/ satisfies the conditions 

\\F{xo)\\y<V, \\F'ixor'\\LiY,x)<(T-\ 



- F'(x2)||L(x,y) < -^^Iki - X2IIX for \\xj - xo\\x < 2r]a~^ , 
and 2L(r^^ri < 1, 

then there exists x & X such that F{x) = and \\x — xo\\x ^ 2r]a^^. 



Proof. The result follows, for example, by applying Theorem 2.1 in 22 with the choices 
R = 2ria^^, uj{xq,R) = LR and u{xo,R) = \LR^. Similar results can be obtained by 
tracking the constants in most proofs of the inverse function theorem, and assuming local 
Lipschitz continuity of F'. □ 

For our purposes, we set X = U^'°°, Y = U~^'°°, F{u) = S'^^^{yqce,F + u), and xq = 0. As- 
suming that 61, 62 are sufficiently small, our previous analysis gives the residual and stability 
estimates 

V = 0f(^Ci5i(52 + |c2(5^) and a = \(j)"p, 

and, in particular, 

2r/(T"^ = 4ci5i52 + 02^1. 



To ensure that i?2^i.oo(0, 2ria^^) remains within the region of differentiability of F, that is, 
to ensure that {yqce,F + u)'t > for ||m'||£cx) < 2ria~^, it is clearly enough to assume that Si 
and 62 are sufficiently small. 

A modification of (35) then allows the choice L = 2(j)pC2 for the local Lipschitz constant. 

Thus, the condition ensuring the existence of a solution yqcc,F of (11) becomes 

ALa-^ri = 64ciC25i52 + 16c^5i < 1, 

which is satisfied, once again, if we assume that and 62 are sufficiently small. An applica- 
tion of Lemma [7] concludes the proof of Lemma [5] 
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